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ABSTRACT 

Using a sample of over 25000 spectroscopically confirmed quasars from the Sloan Digital Sky 
Survey, we show how quasar variability in the rest frame optical/UV regime depends upon rest 
frame time lag, luminosity, rest wavelength, redshift, the presence of radio and X-ray emission, 
and the presence of broad absorption line systems. Imaging photometry is compared with three- 
band spectrophotometry obtained at later epochs spanning time lags up to about two years. The 
large sample size and wide range of parameter values allow the dependence of variability to be 
isolated as a function of many independent parameters. The time dependence of variability (the 
structure function) is well-fit by a single power law with an index 7 = 0.246 ±0.008, on timescales 
from days to years. There is an anti-correlation of variability amplitude with rest wavelength - e.g. 
quasars are about twice as variable at 1000A as 6000A - and quasars are systematically bluer when 
brighter at all redshifts. There is a strong anti-correlation of variability with quasar luminosity 
- variability amplitude decreases by a factor of about four when luminosity increases by a factor 
of 100. There is also a significant positive correlation of variability amplitude with redshift, 
indicating evolution of the quasar population or the variability mechanism. We parameterize 
all of these relationships. Quasars with RASS X-ray detections are significantly more variable 
(at optical/UV wavelengths) than those without, and radio loud quasars are marginally more 
variable than their radio weak counterparts. We find no significant difference in the variability 
of quasars with and without broad absorption line troughs. Currently, no models of quasar 
variability address more than a few of these relationships. Models involving multiple discrete 
events or gravitational microlensing are unlikely by themselves to account for the data. So-called 
accretion disk instability models are promising, but more quantitative predictions are needed. 

Subject headings: galaxies: active - quasars: general - techniques: photometric 
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1. Introduction 

The luminosities of quasars and other active 
galactic nuclei (AGNs) have been observed to 
vary from X-ray to radio wavelengths, and on 
time scales from several hours to many years. 
The majority of quasars exhibit continuum vari- 
ability on the order of 10% on timescales of 
months to years. A minority of AGNs, broadly 
classified as blazars, vary much more dramati- 
cally on much shorter timescales. The mecha- 
nisms behind quasar variability are not known, al- 
though in principle variability is a powerful means 
of constraining models for the energy source of 
AGNs. The most promising models (for non- 
blazar variability) include accretion disk insta- 
bilities (e.g. Rees 1984; Kawaguchi, Mincshigc, 
Umemura, & Turner 1998), so-called Poissonian 
processes such as multiple supernovae (e.g. Ter- 
levich, Tenorio-Tagle, Franco, & Melnick 1992) 
or star collisions (Courvoisier, Paltani, & Walter 
1996; Torricelli-Ciamponi, Foellmi, Courvoisier, & 
Paltani 2000), and gravitational microlensing (e.g. 
Hawkins 1993). Only recently have the various 
models become quantitative enough for meaning- 
ful comparison with observations. A consensus on 
the observational trends with variability is emerg- 
ing, but disagreements remain and even the most 
fundamental relationships need better characteri- 
zation. 

Several dozen studies of quasar optical broad- 
band variability have appeared in the literature. A 
number of the more important studies are summa- 
rized in tabular form by Helfand et al. (2001) and 
Giveon et al. (1999). Most ensemble studies have 
focused on establishing correlations between vari- 
ability (defined in various ways as a measure of the 
source brightness change) and a number of param- 
eters, most importantly time lag, quasar luminos- 
ity, rest frame wavelength, and redshift. Charac- 
teristic timescales of variability range from months 
to years (e.g. Collier & Peterson 2001; Cristiani 
et al. 1996; di Clemente et al. 1996; Smith & 
Nair 1995; Hook, McMahon, Boyle, & Irwin 1994; 
Trcvese et al. 1994). The amplitude of variability 
rises quickly on those timescales, but may slow or 
even level off on longer timescales. 

An anti-correlation between quasar variabil- 
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ity and luminosity was reported by Angione & 
Smith (1972), and confirmed in numerous sub- 
sequent studies (Uomoto, Wills, & Wills 1976; 
Pica & Smith 1983; Lloyd 1984; O'brien, Gond- 
halekar, & Wilson 1988; Hook, McMahon, Boyle, 
& Irwin 1994; Trevese et al. 1994; Cid Fernan- 
des, Aretxaga, & Terlevich 1996; Cristiani et al. 
1996; Cristiani, Trentini, La Franca, & Andreani 
1997; Paltani & Courvoisier 1997; Giveon et al. 
1999; Garcia, Sodre, Jablonski, & Terlevich 1999; 
Hawkins 2000; Webb & Malkan 2000). Such an 
anti-correlation is expected in Poissonian models, 
although complex versions are necessary to explain 
the diversity of the relationship among quasars 
(Cid Fernandes, Sodre, & Vieira da Silva 2000). 

There is strong evidence from multiwavelength 
observations of quasars that variability increases 
with decreasing rest wavelength, which holds over 
a wavelength range spanning at least the ultravi- 
olet to near infrared (Cutri, Wisniewski, Rieke, & 
Lebofsky 1985; Neugebauer, Soifer, Matthews, & 
Elias 1989; Kinney, Bohlin, Blades, & York 1991; 
Paltani & Courvoisier 1994; di Clemente et al. 
1996; Cristiani, Trentini, La Franca, & Andreani 
1997; Giveon et al. 1999; Cid Fernandes, Sodre, & 
Vieira da Silva 2000; Helfand et al. 2001; Trevese 
& Vagnetti 2002). The wavelength dependence 
is related to the observed tendency for quasar 
spectra to become harder (bluer) in bright phases 
(Cutri, Wisniewski, Rieke, & Lebofsky 1985; Edel- 
son, Pike, & Krolik 1990; Giveon et al. 1999; Cid 
Fernandes, Sodre, & Vieira da Silva 2000; Trevese, 
Kron, & Bunonc 2001). The chromatic nature 
of quasar variability is often taken as evidence 
against gravitational microlensing as the primary 
cause of variability (e.g. Cristiani, Trentini, La 
Franca, & Andreani (1997); except see Hawkins 
(1996)), although this may be accounted for if re- 
gions closer to the center are both brighter and 
bluer. 

A correlation of variability with redshift is often 
reported (Cristiani, Vio, & Andreani 1990; Gial- 
longo, Trevese, & Vagnetti 1991; Hook, McMahon, 
Boyle, & Irwin 1994; Trevese et al. 1994; Cid Fer- 
nandes, Aretxaga, & Terlevich 1996; Cristiani et 
al. 1996; Trcvese & Vagnetti 2002) if wavelength 
and luminosity dependencies are not taken into ac- 
count. For a fixed observer timescale, the increase 
of variability with increasing redshift would con- 
tradict the expected 1 + z effect of time dilation. 
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However, it has been shown by Giallongo, Trevese, 
& Vagnetti (1991), Cristiani et al. (1996), and Cid 
Fernandes, Aretxaga, & Terlevich (1996) that the 
inverse wavelength dependence can easily account 
for the uncorrected redshift correlation, since for 
a fixed passband in the observer frame, quasars 
with higher redshifts are detected at shorter wave- 
lengths, which systematically vary at a greater 
amplitude. It is still not clear whether any red- 
shift dependence remains after accounting for rest 
wavelength and luminosity (which is strongly cor- 
related with redshift in flux-limited quasar sam- 
ples). Some studies which have leverage in both 
redshift and luminosity suggest a weak correla- 
tion of redshift and variability (Hook, McMahon, 
Boyle, & Irwin 1994; Cristiani et al. 1996), but 
others show no such effect (Cimatti, Zamorani, & 
Marano 1993; Paltani & Courvoisier 1994; Netzer 
et al. 1996; Cristiani, Trentini, La Franca, & An- 
dreani 1997; Helfand et al. 2001). 

Variability is sometimes found to be correlated 
with radio loudness (Pica & Smith 1983; Smith & 
Nair 1995; Garcia, Sodre, Jablonski, & Terlevich 
1999; Eggers, Shaffer, & Weistrop 2000; Helfand et 
al. 2001; Enya et al. 2002), the equivalent width 
of the H/3 line (Giveon et al. 1999; Cid Fernandes, 
Sodre, & Vieira da Silva 2000), and the presence of 
broad absorption line troughs (Sirola et al. 1998), 
although the results arc not conclusive. No large 
X-ray detected quasar sample has been systemati- 
cally studied for optical variability, but since most 
blazars are X-ray bright, a greater degree of vari- 
ability may be expected from such a sample (e.g. 
Ulrich, Maraschi, & Urry 1997). 

In this paper we present results on a quasar 
variability program using data from the Sloan Dig- 
ital Sky Survey (SDSS, York et al. 2000). A com- 
plementary variability study by de Vries, Becker, 
& White (2003) presents a comparison of the SDSS 
Early Data Release (Stoughton et al. 2002) imag- 
ing photometry with archival photographic plate 
data. One of the goals of the present work is to 
characterize the spectroscopic calibrations of the 
SDSS, in order to examine the spectroscopic vari- 
ability properties of quasars and other objects ob- 
served in the spectroscopic survey. The present 
work uses the broad band fluxes of the spectra 
convolved with the SDSS filter transmission func- 
tions in direct comparison with the imaging pho- 
tometry. This provides photometric data at two 



epochs in three bands for every spectroscopically 
confirmed quasar in the survey - a sample size 
currently of over 25000 quasars. This is by far 
the largest quasar UV/optical variability study to 
date, and it also includes the largest samples of ra- 
dio selected, X-ray selected, and broad absorption 
line quasars ever examined for variability. Our 
goal is to characterize the ensemble dependence of 
variability on many quasar parameters and types, 
on timescales from weeks up to several years. 

We describe the quasar sample drawn from the 
SDSS in § 2. Ensemble measurements of the vari- 
ability are given in § 3. We disentangle the de- 
pendence of variability upon time lag, luminos- 
ity, wavelength, and redshift in §4, show how 
quasar colors change with variability in § 5, and 
look at variability in various quasar subclasses 
in § 6. The implications of the results are dis- 
cussed in § 7, and we conclude in § 8. Through- 
out the paper we assume a flat, cosmological con- 
stant dominated cosmology with parameter values 
tt A = 0.7, Q M = 0.3, and H a = 65km/s/Mpc. 

2. The Quasar Dataset 

2.1. The Sloan Digital Sky Survey 

The Sloan Digital Sky Survey (SDSS) is a 
project to image 10 4 deg 2 of sky mainly in the 
northern Galactic cap, in five broad photometric 
bands (u,g,r, i,z) to a depth of r <~ 23, and to 
obtain spectra of 10 6 galaxies and 10 5 quasars ob- 
served in the imaging survey (York et al. 2000). 
All observations are made with a dedicated 2.5m 
telescope at Apache Point Observatory in New 
Mexico. Images are taken with a large mosaic 
CCD camera (Gunn et al. 1998) in a drift-scanning 
mode. Absolute astrometry for point sources is 
accurate to better than 100 milliarcseconds (Pier 
et al. 2003). Site photometricity and extinction 
monitoring are carried out simultaneously with 
a dedicated 20in telescope at the observing site 
(Hogg, Finkbeiner, Schlegel, & Gunn 2001). The 
imaging data are reduced and calibrated using the 
PHOTO software pipeline (Lupton et al. 2001). In 
this study we use the point-spread function (PSF) 
magnitudes, which are determined by convolving 
the reduced imaging data with a model of the 
spatial point-spread function. The PSF magni- 
tudes are more stable than aperture magnitudes 
for point sources, since they are less dependent 
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upon seeing variations, and because the PSF back- 
ground noise is less within the survey seeing limit 
(which is 1.7 arcseconds). The SDSS photometric 
system is normalized so that the u,g,r,i,z mag- 
nitudes arc on the AB system (Smith et al. 2002; 
Fukugita et al. 1996; Okc & Gunn 1983). The pho- 
tometric zeropoint calibration is accurate to bet- 
ter than 1% (root-mean-squared) in the g, r, and 
i, bands, and to better than 2% in the u and z 
bands, measured by comparing the photometry of 
objects in scan overlap regions. The SDSS image 
reduction and calibration routines have evolved 
throughout the course of the survey and the imag- 
ing runs have been reprocessed accordingly. Thus 
the object imaging magnitudes deemed "best", 
and which we use in this study, may be slightly dif- 
ferent than those used for the spectroscopic target 
selection, although any differences are insignificant 
to the results of this study. Throughout this paper 
we use magnitudes corrected for Galactic extinc- 
tion according to Schlegel, Finkbciner, & Davis 
(1998). 

Objects are selected for spectroscopic follow-up 
as candidate galaxies (Strauss et al. 2002; Eisen- 
stein et al. 2001), quasars (Richards et al. 2002), 
and stars (Stoughton et al. 2002). The spectro- 
scopic targets are grouped by 3 degree diame- 
ter areas or "tiles" (Blanton et al. 2003). For 
each tile, an aluminum plate is drilled with holes 
corresponding to the sky locations of the targets 
along with holes for blank sky, calibration stars, 
and guide stars. The plates are placed at the fo- 
cal plane of the telescope, and optical fibers run 
from the hole positions to two spectrographs, each 
of which accepts 320 fibers allowing for the si- 
multaneous observation of 640 objects. For each 
plate, approximately 500 galaxies, 50 quasars, and 
50 stars are observed. Spectroscopic observations 
generally occur up to a few months, but occasion- 
ally years, after the corresponding imaging obser- 
vations, depending upon scheduling constraints. 
The spectroscopic data for this study come from 
479 spectroscopic plates observed and processed 
through September 2002; 284 of the plates are part 
of the SDSS First Data Release (DR1, Abazajian 
et al. 2003), publicly available since April 2003. 
Seven of the (291) DR1 plates are not included in 
this study, since the DR1 plate list was not final- 
ized until after the sample for this study had been 
gathered. 



2.2. Quasar Target Selection and Sample 
Definition 

Quasar candidates are selected from the imag- 
ing sample by their non-stellar colors from the 
five-band photometry as well as by matching 
SDSS point sources with FIRST radio sources 
(Stoughton et al. 2002). The selection is similar 
to that described by Richards et al. (2002), but 
the formal implementation of this algorithm was 
imposed after the cutoff date for the DR1 quasar 
sample, and much of the post-DRl data used in 
this study. About two-thirds of the candidates are 
confirmed to be quasars from the spectroscopic 
survey. Ultraviolet excess quasars are targeted to 
a limit of i = 19.1 and higher rcdshift quasars 
are targeted to i = 20.2. These criteria give a 
sample that is estimated to be over 90% complete 
(Richards et al. 2002). Additional quasars are tar- 
geted as part of the SERENDIPITY and ROSAT 
classes (Stoughton et al. 2002) or (incorrectly) as 
stars. 

Quasars are identified from their spectra us- 
ing a combination of both automated classifica- 
tion (about 94%) and manual inspection (about 
6%) of those objects flagged by the spectroscopic 
pipeline as being less reliably identified. For the 
purposes of this study, we define "quasar" to mean 
any extragalactic object with broad emission lines 
(pipeline measured full width at half maximum ve- 
locity width of > lOOOkms -1 ) regardless of lumi- 
nosity. The definition thus includes objects which 
are often classified as less luminous types of ac- 
tive galactic nuclei (AGNs) rather than quasars, 
and excludes AGNs without strong broad emis- 
sion lines such as BL Lacs and some extreme broad 
absorption line quasars. To assemble our sample, 
we extract relevant data for all point sources from 
the SDSS database. Data from only one imaging 
and one spectroscopic epoch are used per object 
to avoid giving extra weight to any object. Ex- 
tended sources are rejected because they compli- 
cate the spectrophotometric recalibration (§2.4), 
and their spectra are likely to be seriously con- 
taminated by host galaxy light. Objects with bad 
spectra (defined to be those with significantly long 
unprocessable portions of spectrum) are rejected. 
Those remaining objects spectroscopically identi- 
fied as stars are used to refine the spectroscopic 
calibration (see §2.4). The remaining 25710 ob- 
jects identified to be quasars become part of the 
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quasar variability sample. 

A catalog of quasars found in the SDSS DR1 
dataset is described by Schneider et al. (2003). Of 
the 16713 listed DR1 quasars, 14705 are included 
in this study. The difference can be attributed to 
two factors. First, the definitions of "quasar" are 
slightly different - Schneider et al. (2003) impose 
an absolute magnitude limit of Mi < —22.0, and 
they make no distinction between point-like and 
extended sources. Second, data from only 284 of 
the 291 DR1 spectroscopic plates are included in 
this study. Over half (57%) of the quasars in this 
study are contained in the DR1 sample, and we ex- 
pect that our results on quasar variability would 
be similar (albeit more noisy) for the full DR1 
sample, except possibly for the extended sources. 

Absolute magnitudes in the rest frame i band 
Mi, are calculated for each quasar using dered- 
dened observed i band PSF magnitudes and as- 
suming a power law spectral energy distribution 
fx cx A Qa , with a wavelength slope of ax = —1.5. 
Detailed K-corrections are not yet available for 
quasars in the SDSS photometric system. How- 
ever, estimated K-corrections using the compos- 
ite spectrum from Vanden Berk et al. (2001) are 
consistent with the simple power law assumption, 
and the differences are usually no greater than 0.1 
magnitude at any rcdshift. The differences are not 
significant here since the data are coarsely binned 
(by statistical necessity) when we examine abso- 
lute magnitude trends. 

2.3. SDSS Spectroscopy and Its Calibra- 
tion 

Spectra are obtained in three to four consec- 
utive 15-minute observations. There are 32 sky 
fibers, 8 spectrophotometric standard stars and 8 
reddening standards stars observed on each plate 
to help with calibration of the remaining 592 sci- 
ence targets. Spectral reductions and calibrations 
are done using the SDSS SPECTRO 2D pipeline 
(Stoughton et al. 2002). The 8 spectrophotomet- 
ric calibration stars are chosen to approximate the 
standard F0 subdwarf star BD + 17°4708, and are 
used by the 2D pipeline for absolute spectral flux 
calibration and dereddening due to Galactic ex- 
tinction. The 2D pipeline also calculates synthetic 
spectroscopic magnitudes by convolving the cali- 
brated spectra with SDSS g, r, and i filter trans- 
mission curves, assuming 1.2 airmasses of extinc- 



tion (the spectra do not cover the entire wave- 
length ranges of the u and z bands). By using 
these synthetic spectral magnitudes, we obtain a 
second photometric data point for every spectro- 
scopically observed quasar. Additionally, a signal- 
to-noise (S/N) parameter is calculated for each 
of the three bands by convolving the spectral er- 
ror with the transmission curves and dividing that 
into the corresponding convolved flux. The spec- 
tral magnitude signal-to-noise is essential for char- 
acterizing magnitude difference uncertainties (see 
§2.4). 

2.4. Refinement of Spectroscopic Calibra- 
tion 

The differences between the spectroscopic and 
imaging PSF magnitudes Am = m s — m p for all 
of the spectroscopically confirmed stars in each of 
the g, r, and i bands are shown in Fig. 1. Well 
calibrated data should center around zero mag- 
nitude difference, and there should be no trend 
with magnitude, except for larger uncertainties 
at fainter magnitudes. There are clearly system- 
atic differences in the magnitudes derived from the 
spectra and from the images. While the initial 
spectroscopic calibration is more than adequate 
for the primary purposes of the SDSS, namely ob- 
ject identification and reliable redshift measure- 
ment, variability studies require more careful cal- 
ibrations. Fortunately, the magnitude difference 
trends with PSF magnitude can be easily under- 
stood and are almost entirely correctable. 

There are three primary sources for the mag- 
nitude difference discrepancy: the inclusion of ob- 
jects with bad PSF magnitudes, an aperture ef- 
fect relating the finite fiber diameters to the PSF 
magnitudes, and what may be a very small but 
significant sky under-subtraction in the spectro- 
scopic data. Occasionally, point sources in the im- 
ages can have poor photometry if they are closely 
blended with other objects, occur where the seeing 
has changed very rapidly, or lie where there may be 
other problems in the imaging data. The long tails 
in the histograms of Fig. 1 are populated mainly 
by the measurements of these objects. Because the 
objects will have unusual measured colors they are 
sometimes selected as high-rcdshift quasar candi- 
dates for spectroscopic follow-up, which turn out 
to be normal stars upon examination of the spec- 
tra. Therefore, spectroscopically confirmed stars 
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which were selected as high-redshift quasar candi- 
dates are removed from the stellar data set for 
spectrophotometric refinement. Late- type stars 
identified by the spectroscopic pipeline are also re- 
jected because they are often variable. 

The median Am offset from zero is simply an 
aperture effect wherein the 3" spectroscopic fibers 
subtend a smaller fraction of the total object im- 
age than the point-spread function used to mea- 
sure the PSF magnitudes in the imaging data. The 
spectroscopic fiber flux density to PSF flux density 
ratio is nearly a constant (but somewhat depen- 
dent upon seeing at the spectroscopic epoch, see 
below), so the magnitude difference will also be a 
nearly (non-zero) constant. 

The downward trend of Am with PSF magni- 
tude seen for each band is most easily accounted 
for by a small overestimation of the flux density 
in the spectroscopic data, possibly caused by a 
slight under-subtraction of the sky level. Tests 
of the imaging photometric calibration show that 
the effect is not likely to be caused by sky over- 
subtraction in the imaging data. Further tests will 
have to be done to determine the cause with cer- 
tainty. The correction for a flux density overesti- 
mation combined with a fiber aperture correction 
give 

Am = m s -m p = -2.51og(/ s // p ) 
= -2.5 log(a(/ p + &)//„) 
= -2.51og(a + a610 (m ^ c ' )/2 - 5 ) , (1) 

where a is the aperture correction, b is the cor- 
rection for the flux density offset, and C is the 
zeropoint constant used in converting flux density 
to magnitude. Assuming a and b arc constants, 
the function has two adjustable parameters, and 
fits to the data provide reasonably good descrip- 
tions of the Am vs. m p trend. However, in or- 
der to account for any other effects, expected or 
unexpected, we use the following more flexible 3 
parameter function 

Am = Arno — exp((m p — mo)/m*) , (2) 

where Amo, m , and m* are constants to be deter- 
mined from the fits to the data. For example, the 
Malmquist bias (Malmquist 1924) will add to the 
magnitude difference approximately as the square 
of the PSF magnitude uncertainty, <j p 2 . Since a p 
ranges from about 0.01 to 0.05, the Malmquist bias 



is expected to affect the magnitude difference by 
at most a few percent. In the absence of this or 
other higher order effects, the 3-parameter func- 
tion would almost exactly reproduce the 2 param- 
eter logarithmic fit. The 3-parameter functions 
are fit to the data in each band separately, then 
subtracted from the magnitude differences. 

After this correction, offsets from zero remain 
for the mean Am values for stars on the same 
plate, in excess of those expected from statistical 
uncertainty. These plate-to-plate offsets arc due 
to differences in the spectral energy distributions 
of the stars used as initial spectrophotometric cal- 
ibrators relative to BD + 17°4708 and due to dif- 
ferences in the seeing at the epochs of the plate 
observations. The former effect applies to "half 
plates" corresponding to the two sets of 320 fibers 
running separately to each spectrograph. To cor- 
rect for these effects, we reject stars outside of the 
99% confidence envelope resulting from the Am 
corrections described above, and work only with 
half plates which have at least 5 remaining stars 
(the average number is about 20). The median 
Am offsets are calculated for each half plate and 
subtracted from the values for each of the stars 
observed with that half plate. 

The final corrected Am distributions are shown 
in Fig. 2 as a function of spectral S/N. The width 
of the stellar Am distribution is correlated with 
magnitude, but the better correlation is with spec- 
tral S/N. The reason for this is that while magni- 
tude and S/N are correlated, it is the S/N which is 
directly related to the quality of a spectrum. As 
a whole, the 68.3% confidence half-width (nomi- 
nally 1 standard deviation) is « 0.08 in each band 
at a spectroscopic S/N of 10, which is a substan- 
tial improvement over the initial widths of « 0.13. 
Fits to the 68.3% confidence half-width as a func- 
tion of spectral S/N are shown in Fig. 2, for which 
we used a function of the form 

&s/n = ao + aiexp(a 2 5/A) , (3) 

where ao, a\, and a 2 are constants. These fits are 
used as statistical measurement uncertainties for 
the quasars (see § 3). 

The same spectrophotometric corrections ap- 
plied to the stars are also applied to the quasars. 
The resulting distribution of quasar magnitude 
differences as a function of spectral S/N, and the 
histograms of magnitude differences are shown in 
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Fig. 3. The mean corrected magnitude differences 
for the quasars are 0.002, -0.004, and -0.011 for 
the g, r, and i bands respectively. These values 
are small compared with the measurement uncer- 
tainties derived from the stars. It is possible that 
small differences in the spectral energy distribu- 
tions (SEDs) of stars and quasars affect the re- 
calibration of the quasar photometry. However, 
any effect is expected to be quite small since the 
same filter transmission curves are used for both 
the imaging and spectroscopic photometry, and 
the majority of the stars used for the recalibration 
were selected as quasar candidates in the color- 
selected survey, which guarantees that the SEDs 
are very similar. The 68.3% confidence limit half- 
widths of the g, r, and i magnitude difference dis- 
tributions are 0.134,0.119, and 0.114 respectively 
at a S/N of 10, substantially larger than those 
of the stars. The stars and quasars were selected 
to be point sources, observed simultaneously with 
the same instrument, and often were selected with 
the same algorithm. The larger magnitude differ- 
ences among the quasars, therefore, demonstrate 
the variable nature of the quasars in the sample. 
The following sections quantify the variability and 
its dependence upon many quasar parameters. 

3. Ensemble Variability Measurement 
The Structure Function 

The magnitude difference histograms from the 
previous section show that the quasars are signif- 
icantly more variable as a class than the stars. 
Assuming no stellar variability, we can use the 
distribution of the stellar magnitude differences 
to quantify the statistical measurement uncertain- 
ties. Removing the width of the stellar magnitude 
difference distribution in quadrature, the average 
quasar magnitude differences (at a spectral S/N 
of 10) due to variability in the sample are 0.103, 
0.086, and 0.080 in the g, r, and i bands respec- 
tively Measurement uncertainties must be taken 
into account because they are comparable to the 
values of the variability itself. The large sizes of 
the samples (both the quasars and the compari- 
son stars) allow the measurement uncertainty to 
be effectively removed. 

We first show the absolute values of the mea- 
sured quasar magnitude variations, uncorrected 
for measurement uncertainty, in Figures 4 through 



7, as a function of quasar rest frame time lag 
(Fig. 4) , absolute magnitude in the rest frame i 
band (Fig. 5), rest frame wavelength (Fig. 6), and 
redshift (Fig. 7). Data in each of the three photo- 
metric bands are shown separately. Average val- 
ues in a set of bins are also shown. Because flux 
densities in the Ly a forest region are not repre- 
sentative of the true quasar flux, we have omitted 
data in each band at redshifts beyond which the 
Lya forest covers the band: z = 2.5, 4.75, and 
6.0 for the g, r, and i bands respectively. The 
number of measurements rejected for each band 
due to the Lya forest are 742, 45, and for the 
g, r, and i band respectively. The figures show 
that there are several apparent correlations even 
in the uncorrected data. In particular, the aver- 
age magnitude difference increases with time lag 
in all three bands, and the magnitude difference 
decreases with more negative absolute magnitude 
(decreases with luminosity). No trends are ap- 
parent at this stage between magnitude difference 
and rest wavelength or redshift. Again, it is im- 
portant to account for measurement uncertainty 
before making any claims about the dependence 
of variability on any parameter. 

The definition of variability used here is a sta- 
tistical measure of the magnitude difference, tak- 
ing into account measurement uncertainty. The 
first application is to the dependence of variabil- 
ity on rest frame time lag - the so-called structure 
function. Historically, the structure function has 
been the primary measure of variability for stud- 
ies of both individual quasars and quasar ensem- 
bles. For individual quasars with multiple sam- 
pling epochs, the structure function is comprised 
of the values of the magnitude differences for each 
pair of time lags in the data set, and it is closely re- 
lated to the autocorrelation function (e.g. Simon- 
ctti, Cordes, & Heeschen 1985). In the ensem- 
ble case, here with only two sampling epochs, the 
structure function is simply the average value of 
the magnitude difference for all objects with the 
same (or nearly the same) time lags. The error 
analysis is simpler in the ensemble case since all 
of the data points are independent. 

We define the ensemble variability, V, of a set 
of quasars as 

V = (|(|Am|) 2 - (a 2 s/N )f , (4) 
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where Am is the measured magnitude difference, 
(Js/n is the statistical measurement uncertainty of 
Am (as a function of spectral S/N) derived from 
the fits to the star measurements in §2.4, and 
the brackets denote average quantities. The av- 
erage absolute value of the magnitude difference, 
along with the scaling factor of 7r/2, is more robust 
against the presence of outliers in the data than 
the average of the square of the differences. The 
values V function of rest frame time lag At, 
define the structure function, V(At). The same 
relation has been used for the structure function 
in previous variability studies (e.g. di Clemente et 
al. 1996). 

The binned structure function for all of the 
quasars in the sample for each of the three photo- 
metric bands is shown in Fig. 8 with logarithmic 
axes. The bins were chosen to have equal intervals 
in logarithmic rest frame time lag, and to have 
reasonably large numbers of objects. The number 
of objects per bin range from 241 for the shortest 
time lag bin covering 7 to 11 days in the g band, to 
7919 for the i band bin covering time lags from 111 
to 176 days. Quasars with magnitude differences 
larger than 0.75 - just over 5 times the la width of 
the distribution for quasars - were rejected from 
the analysis in order to remove outliers. This step 
removes about 1% of the quasars, which is more 
than would be expected for a truly normal dis- 
tribution. In a related paper (Ivezic et al. 2003), 
the distribution of ultraviolet-excess quasar mag- 
nitude differences, including very large differences, 
is discussed in more detail. The apparently highly 
variable quasars may be optically violent variables 
and are valuable for follow-up studies, but the fo- 
cus here is on "typical" quasar variability. The 
error bars were determined by propagating the 
root-mean-square errors a 2 , in the average mag- 
nitude difference and measurement uncertainty in 
quadrature 

a(V) = iy- 1 v /^ { | Am | > 2 CT2({ | Am | ))+<72({(T 2 /jv>) . 

Two trends are obvious: first the structure func- 
tion increases as a function of time lag - the mag- 
nitude differences are greater the longer the time 
between measurements. Second, the amplitude of 
variability is greater in the g (bluest) band than in 
the others, and the r band amplitude is generally 
greater than in the i (reddest) band. This is the 
variability anticorrelation with wavelength found 



in a several previous studies. We will quantify 
the wavelength dependence explicitly in § 4.4, ac- 
counting for the dependence on other parameters. 
For the purposes of the remainder of this section, 
the clear wavelength dependence means that the 
analysis of the structure function will address the 
three bands individually. 

The correlation of variability with time lag has 
been found in numerous previous studies, however 
the form of the structure function has remained a 
topic of debate. We fit the binned structure func- 
tions with the two most common parameteriza- 
tions. The first is a power law (e.g. Hook, McMa- 
hon, Boyle, & Irwin 1994) 

"(Ar)-(^)' , (6) 

where Ar and 7 are constant parameters to be 
determined. This will appear as a straight line 
in a log-log plot such as Fig. 8. The second is an 
asymptotic function - a constant minus an expo- 
nential - which is the most common parameteri- 
zation of the structure function (e.g. Bonoli et al. 
1979; Trevese et al. 1994; Hook, McMahon, Boyle, 
& Irwin 1994; Enya et al. 2002) 

V{At) = V (l - e^) , (7) 

where again Vb and Atq are to be determined. 
The "time scales", Aro, whatever their values, are 
simply parameters of the functions to be fit to the 
data, and can not necessarily be directly compared 
with physical characteristic time scales, such as 
those associated with accretion disks, star bursts, 
or gravitational lens dynamics. 

Parameter values, uncertainties, and x 2 values 
for each of the functions in each of the photo- 
metric bands are given in Table 1. Based on the 
X 2 values, the functional form that best fits the 
structure function in each band is a power law. 
The power law slopes in each of the three bands 
( 5 )- 0.293 ± 0.030,0.336 ± 0.033, and 0.303 ± 0.035 
for g, r, and i respectively - are consistent with 
each other within one standard deviation of the 
difference. The power law scale factors (where the 
structure functions would have a value of one) are 
not well constrained, mainly because the observed 
time lag only extends to about 700 days. The 
shape of the structure function at much longer 
time lags is sometimes observed to "flatten" some- 
what (Cristiani, Trentini, La Franca, & Andreani 
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(1997); Hook, McMahon, Boyle, & Irwin (1994); 
Trevese et al. (1994); see Hawkins (2002) for a 
counter-example), but the data at long time lags 
do not yet favor one parameterization over an- 
other. In any case, what we can say from this 
study is that a 2-parameter power law is a good 
description of the data - and a better description 
than a 2-parameter exponential - up to time lags 
of about 2 years. As the SDSS proceeds, the range 
of time lags will eventually reach up to about 5 
years. The rest frame time sampling will continue 
to improve on all scales, and the power law form 
can be even more stringently tested on longer time 
scales. 

The wavelength dependence of the structure 
function becomes clearer at longer time lags (and 
should become even clearer by survey end), as 
would be expected from either a power law or ex- 
ponential fit. The distributions of quasar proper- 
ties - e.g. luminosity and rest frame wavelength 
- vary with rest frame time lag because of survey 
selection effects and the dependence of these prop- 
erties on redshift. In the next section we disentan- 
gle the dependence of variability on four primary 
quasar parameters. 

4. Variability Dependence on Time Lag, 
Luminosity, Wavelength, and Redshift 

4.1. Selection Function 

The structure function calculated in the pre- 
vious section describes the variability of the full 
data set with respect to time lag. However, vari- 
ability is almost certainly also a function of quasar 
luminosity and rest frame wavelength, and possi- 
bly redshift. In order to separate the dependence 
of variability on multiple parameters, the selec- 
tion biases must be taken into account. Even the 
structure function may not give the true depen- 
dence of variability on rest frame time lag. The 
set of quasars within a narrow range of time lags 
will be populated with objects with wide ranges 
of the other parameters. In addition, since rest 
frame time lag is dependent upon redshift, as are 
the other parameters, the distributions of quasar 
parameter values will be correlated with the time 
lag. For example, high redshift quasars will gener- 
ally have shorter time lags in the rest frame than 
lower redshift quasars. 

The selection function - that is, the region of 



parameter space occupied by the data set - is 
shown in projected planes in Fig. 9. A number 
of artificial correlations are evident and are due 
both to the survey selection criteria, and to the 
dependence of luminosity, rest frame time lag, and 
rest wavelength on redshift. Variability informa- 
tion can obviously only be obtained in the regions 
of the four dimensional parameter space contain- 
ing a statistically sufficient number of objects. 

In order to determine how variability is related 
to a single parameter, the space was divided into 
small regions in three dimensions, and the vari- 
ability calculated as a function of the remaining 
parameter in each of the slices. The condition 
that there be enough quasars to reliably measure 
variability was the primary limiting factor for the 
bin sizes. For each parameter there is then a set 
of variability relations, each set representing the 
results of restricting the ranges of the other 3- 
parameters. As shown in the remainder of this 
section, in most cases variability trends are clear 
even in independent restricted data sets. 

If wc make the assumption that the equa- 
tions describing the multi-parameter dependence 
of variability are separable, the results from each 
of the slices may be scaled in the single parameter 
ranges where they overlap, in order to find the 
variability dependence upon a single parameter. 
That is, the form of the variability dependence is 
assumed to be 

V(At, A, z) = v(At) x v(Mi) x v(X) x v(z) , (8) 

where V is calculated as in the previous sec- 
tion. While this form is not necessarily correct, it 
greatly simplifies the analysis, and the relatively 
simple relations found for each parameter suggest 
that it is not far from reality. In the following 
subsections, we show the unsealed variability rela- 
tions for single parameters in the independent re- 
stricted data sets, then show the results after scal- 
ing the independent sets together assuming equa- 
tion 8. The scaled relations are fit with relatively 
simple descriptive functions for each parameter. 

4.2. Time Lag 

We focus first on the dependence of variability 
on rest frame time lag, independent of the other 
3 parameters. The full quasar sample was first 
separated into 6 redshift bins, each with an equal 
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spacing in logarithmic 1 + z. The redshift bin 
sides are: 0.185, 0.499, 0.895, 1.395, 2.028, 2.829, 
and 3.840. The quasar sample in each redshift bin 
was then divided into two halves separated at the 
median absolute magnitude of the quasars in the 
redshift bin: M hmedian = -22.96, -24.07, -25.24, 
-26.07, -26.73, and -27.09. Taking each photo- 
metric band separately for the quasars in each of 
the 12 redshift/absolute magnitude bins restricted 
the quasar rest wavelengths to small ranges. The 
procedure produced 36 independent data sets con- 
fined to small ranges of redshift, absolute magni- 
tude, and rest wavelength, but unrestricted with 
respect to rest frame time lag. The variability am- 
plitude and uncertainty as a function of time lag 
were determined as in § 3, for the quasars in each 
of the 36 data sets independently. The time lag 
bins were set to have a constant logarithmic time 
width, as in § 3, but with twice the width to ac- 
commodate the smaller number of objects per bin. 
The rest frame time lag bin sides, in days, are: 
7.0, 17.6, 44.3, 111.1, 278.7, and 699.2. Unphysical 
(imaginary) values of the variability amplitude oc- 
curred in a small number (3) of cases in which the 
number of quasars was relatively small. In most 
cases when this occurs the number of quasars is 5 
or fewer. For all further analysis, binned data sets 
containing fewer than 10 quasars or which produce 
imaginary values of the variability amplitude are 
rejected. For each of the 36 data sets, the vari- 
ability with respect to rest frame time lag is an 
independent structure function over which the ab- 
solute magnitude, rest wavelength, and redshift do 
not vary greatly. The results are shown in Fig. 10. 
Results for the three photometric bands, corre- 
sponding to restricted ranges of rest wavelength, 
are shown in separate panels. The average red- 
shift of the quasars contributing to each structure 
function is indicated by color, with redder colors 
representing progressively higher redshifts. The 
structure functions containing the more luminous 
halves of the quasar data sets are shown with solid 
points, and the less luminous halves with open 
points. 

It is clear from Fig. 10 that variability is an in- 
creasing function of rest frame time lag at all red- 
shifts, absolute magnitudes, and rest wavelengths. 
Two other trends can also be seen: the less lu- 
minous quasars vary more than their more lumi- 
nous counterparts (nearly all of the open points 



lie above the solid points within the same red- 
shift bin); and quasars vary more at shorter wave- 
lengths, confirming what was shown by the unre- 
stricted structure functions in § 3. 

Under the assumption that the variability as a 
function of time lag is separable from the other de- 
pendencies, the individual structure functions can 
be scaled together in the time lag regions where 
they overlap. Using as a reference a structure 
function near the middle of the redshift, luminos- 
ity, and wavelength distributions, all 36 structure 
functions (excluding bins with too few objects) 
were scaled such that the sum of the products of 
the amplitudes and the time lag bin widths (the 
areas under the curves) were equal. The scaled 
points are shown in the last panel of Fig. 10 along 
with the best fit power law. The parameter fits 
to a power law and exponential are given in Ta- 
ble 2. A 2-parameter power law provides a very 
good fit and is better than the asymptotic (ex- 
ponential) form. The power law fit has a slope 
of 7 = 0.246 ± 0.008, which is comparable to 
but shallower than the values found for the un- 
restricted structure functions in §3. Scaling the 
data points will change the characteristic time 
scale of the function (the time lag at which the 
power law amplitude would be unity), but not 
the power law slope. That the slope is relatively 
close to those found for the unrestricted structure 
functions is due to the offsetting variability depen- 
dencies on luminosity and rest wavelength. From 
Fig. 9, longer time lags are statistically populated 
by more luminous objects, which vary less, but 
at shorter wavelengths, where the variability is 
greater. The significance of the power law slope 
in relation to variability models is discussed in § 7. 

4.3. Absolute Magnitude 

The luminosity (absolute magnitude) depen- 
dence of variability is separated from the other 
parameters in a similar way to the time lag in the 
previous subsection. The full quasar sample was 
separated into the same 6 redshift bins, and each 
separate quasar sample was further divided into 
two halves at the median rest frame time lag of the 
quasars in the redshift bin. The median rest frame 
time lags, in days, for each redshift bin are: 249.9, 
192.0, 157.1, 122.9, 100.5, and 73.6. Again con- 
sidering the three photometric bands separately 
restricted the rest wavelengths of each data set 
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to small ranges. This produced 36 independent 
quasar subsamples unrestricted with respect to ab- 
solute magnitude. The variability amplitude as a 
function of absolute magnitude, in bins one magni- 
tude in width, is shown for each quasar subsamplc 
in Fig. 1 1 . The variability amplitude is an increas- 
ing function of absolute magnitude (brighter ob- 
jects vary less) for nearly every subsample. Also 
seen in Fig. 11 are the time lag dependence - vari- 
ability amplitudes are greater at longer time lags 
(filled symbols) than short time lags (open sym- 
bols) - and a wavelength dependence, seen most 
easily in the amplitude differences between the g 
and i bands. This is a clear demonstration that 
the well-known luminosity-variability anticorrela- 
tion is not simply due to time lag or rest wave- 
length selection effects. 

The data sets were again scaled such that that 
the areas under the curves were equal. In each 
redshift bin, the six individual sets (3 wavelength 
and 2 time lag bins) were scaled to have identical 
sums of the product of the absolute magnitude bin 
width and variability amplitude. Then proceed- 
ing from the lowest to highest rcdshifts, the scaled 
redshift sets were rescaled to the adjacent redshift 
set so that the areas under the curves were equal 
in regions where the absolute magnitude coverage 
overlapped. The scaled data points are shown in 
the last panel of Fig. 11. A straight line can be fit 
to the data points, but such a description is un- 
physical since at large luminosities (large negative 
absolute magnitudes) the function becomes neg- 
ative. In so called Poissonian or discrete-events 
models the relative luminosity variability, SL/L is 
expected to vary with luminosity as SL/L oc L~@, 
where (3 = \ in general (e.g. Cid Fernandes, Sodre, 
& Vieira da Silva 2000). This relationship trans- 
lates into the absolute magnitude form 

v{Mi) oc lO^/ 2 ' 5 . (9) 

This function and one in which (3 is held fixed at 
0.5 were fit to the data. Both fits are shown in 
Fig. 11. The logarithmic equation, with a best fit 
of f3 = 0.246 ± 0.005, fits the data as well as a 
straight line, and avoids the problem of negative 
values. The Poissonian prediction of (3 = 0.5 gives 
a poor fit and is clearly inconsistent with the data. 
Scaling the individual data sets, which accounts 
for arbitrary contributions from variability depen- 
dencies on time lag, rest wavelength, and redshift, 



does not change the value of (3. We will discuss 
the implications for Poissonian models further in 
§7. 

4.4. Rest Wavelength 

The rest wavelength dependence of variability 
was isolated for subsamples selected to cover small 
ranges in redshift, time lag, and absolute magni- 
tude. The redshift and absolute magnitude bins 
are the same ones used to isolate the time lag de- 
pendence (12 separate bins, see §4.2). The data 
in each of these bins were divided into three sepa- 
rate samples in the plane of absolute magnitude 
and time lag, according to the following lines: 
Mi = 0.085Ar- 31.25 and M, = 0.024At - 31.42. 
This gives 36 independent data sets covering small 
ranges of redshift, absolute magnitude, and time 
lag. All but one of the data sets contain enough 
objects to compute reliable variability measure- 
ments. Each data set samples three separate rest 
wavelength points given by the effective rest wave- 
lengths of the three photometric bands. The vari- 
ability amplitude as a function of the rest frame 
wavelength for each set is shown in Fig. 12. The 
average rest frame time lag in each set is color 
coded (longer time lags arc redder), the more lumi- 
nous half of a redshift/time lag bin is shown with 
solid points, and each redshift subset is shown in 
a separate panel. 

In most cases, the variability amplitude de- 
creases with wavelength as expected from previous 
analysis. The cases in which the opposite happens 
occur at short time lags and very low or very high 
redshifts, but there are too few cases to make any 
claims about deviations from the general trend. 
The time lag and luminosity dependencies are also 
evident from the figure. 

The data points were scaled in a manner similar 
to that in §4.3. For each redshift bin, the six sets 
of points (3 time lag bins and 2 absolute magni- 
tude bins) were scaled to the same area under the 
curves and then, moving from low to high redshift, 
the points were rescaled to match the appropriate 
area in the adjacent redshift bin. In this case, 
since the rest wavelength bin limits are not equal 
for the separate redshift bins, the three wavelength 
points in each set were connected by two straight 
lines, and the area under the lines was calculated 
in the regions where they overlapped the wave- 
lengths of the adjacent redshift bin points. The 
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scaled points as a function of rest wavelength are 
shown in Fig. 13. If the variability is due to a sim- 
ple change in the index of a single power law, we 
would expect the wavelength dependence to be 

u(A) = -2.5Aa A log(A) + C , (10) 

where Aa\ is the difference in the wavelength 
power law index, and C is a constant related to 
the pivot wavelength (where the two power laws 
intersect) presumed to be much longer than the 
observed wavelengths. Figure 13 is plotted with 
a logarithmic wavelength axis, and so the relation 
would be seen as a straight line. A single straight 
line is an adequate fit from the shortest wave- 
lengths up to about 4000A, but does not account 
well for the longer wavelength end. Contamina- 
tion from host galaxy light at longer wavelengths 
would have the opposite effect - cause the vari- 
ability to decrease even faster with wavelength. A 
3-parameter exponential function, although phys- 
ically unmotivated, fits the data well 

v(X) = a exp(-A/A ) + Oi , (11) 

with parameter values a>o = 0.616 ± 0.056, Xq = 
988 ± 60A, and ai = 0.164 ± 0.003. 

4.5. Redshift 

The redshift dependence of variability is more 
difficult to isolate from the time lag, absolute mag- 
nitude, and wavelength dependencies. The rea- 
son for this can be seen from inspection of Fig. 9. 
For example, samples restricted to a narrow range 
of rest wavelengths will have 3 independent red- 
shift intervals (taking the three bands separately) , 
but the corresponding absolute magnitude ranges 
for the redshift intervals may not overlap appre- 
ciably. So to isolate the redshift dependence, we 
first found a region of wavelength-absolute mag- 
nitude space which is covered by quasar data in 
all three bands. Fig. 14 shows the wavelength- 
absolute magnitude plane and the selected re- 
gion which is bounded by a triangle with cor- 
ners (A, Mi) : (1250, -27.0), (1250, -29.4), and 
(3491, —23.8). Outside of this region quasar data 
are generally available for only 2 or fewer of the 
bands. The data contained in this region were 
then separated into 16 wavelength bins, shown in 
Fig. 14, at intervals of 100, 200, 400, or 800A, de- 
pending on the number of objects contained in 



the bin. For each bin of wavelength separated 
data, objects were selected from a single range of 
time lags chosen so that the average absolute mag- 
nitudes, time lags, and wavelengths were about 
equal for data in each of the three photometric 
bands. Taking each of the photometric bands sep- 
arately for a restricted data set gives a wide range 
of redshifts, while keeping the ranges of time lag, 
absolute magnitude, and rest wavelength nearly 
constant. 

The variability amplitudes for all 16 data sets 
(each with three redshift points) are shown in 
Fig. 15. The number of objects contributing to 
each point ranges from 65 to 1150, with a mean 
number of 239. Lines connect the points belonging 
to data sets with nearly the same parameter values 
but at different redshifts. The color of the points 
and lines corresponds to the average absolute mag- 
nitude of the data set, with bluer colors represent- 
ing brighter absolute magnitudes. The point sizes 
correspond to the average time lag of the objects 
contributing to the points. The absolute mag- 
nitude correlation, discussed above, is evident in 
Fig. 15, but it is partly counteracted by the nearly 
monotonically increasing rest wavelength with av- 
erage absolute magnitude (seen from Fig. 14), and 
the generally longer average time lags associated 
with fainter average absolute magnitudes. 

What is of interest here is the dependence of the 
variability on redshift. The results are fairly noisy 
and it is difficult to detect any clear trend with red- 
shift. The sets of points were scaled by matching 
the areas under the curves of adjacent data sets, 
as in the previous subsections, starting with the 
sample with the shortest rest wavelengths. The 
scaled data points are shown in Fig. 16. There is 
a correlation between the scaled variability and 
redshift - quasars appear to be more variable at 
higher redshifts. The Spearman rank correlation 
probability that the points are uncorrelated is less 
than 10 -4 , even after accounting for the reduction 
of the number degrees of freedom by the number 
of restricted data sets (16). A straight line fit to 
the data points (linear in redshift and variability 
amplitude) gives 

v(z) = (0.019 ± 0.002)z+ (0.037 ± 0.005) . (12) 

The correlation, although significant, is weak 
enough that it could easily have gone unnoticed in 
previous variability studies, especially since most 
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of them suffer from a lack of sufficient sample over- 
lap needed to test the redshift relationship inde- 
pendently of other parameters. The redshift evolu- 
tion of variability would have serious consequences 
for a number of currently proposed models. If 
the effect is intrinsic, the quasar population or 
the variability mechanism is changing over time. 
External causes are also possible, such as gravita- 
tional microlensing which may increase with red- 
shift since more potential lenses would be avail- 
able. The variability correlation with redshift is 
discussed further in §7. 

5. Color Dependence 

Evidence from previous ensemble studies (Trevese, 
Kron, & Bunonc 2001; Givcon et al. 1999; Edcl- 
son, Pike, & Krolik 1990) suggests that the spec- 
tral energy distribution of quasars becomes harder 
(bluer) in bright phases. Indirect evidence also 
comes from the fact that there is a strong wave- 
length dependence upon variability (§4.4). This 
could happen for example if the index of a power 
law component of the continuum changes with 
luminosity (Trevese & Vagnetti 2002). 

Quasar colors are a strong function of redshift 
(Richards et al. 2001), since various spectral fea- 
tures move into and out of the photometric band- 
passes with redshift changes. A pure power law 
spectrum would have a single set of colors indepen- 
dent of redshift. The observed quasar color struc- 
ture is mainly due to the presence of strong emis- 
sion features, especially broad Fen complexes, as 
well as the Lya forest. Figure 17 shows the av- 
erage imaging photometric colors of quasars as a 
function of redshift in two samples selected to be 
either brighter or fainter by at least 3 standard de- 
viations in at least one of the g, r, i bands relative 
to the spectrophotometry. We use the imaging 
photometric colors rather than the spectroscopic 
since they are more precise, and two more colors 
are available. Although the color differences are 
small (<~ 0.03), at most rcdshifts up to at least 
2.5 and for each color the bright phase sample is 
bluer than the faint phase sample. Also shown 
in Fig. 17 are the color differences of the bright 
phase minus the faint phase samples as a function 
of redshift. Both the binned and average color 
differences are shown. The color differences are 
increasingly larger for shorter wavelength bands, 



i.e. quasars in bright phases are bluer than those 
in faint phases, and they are even bluer at shorter 
wavelengths. 

That the color change persists at high red- 
shift also indicates that it cannot be accounted for 
solely by a non- variable red spectral component, 
such as the quasar host galaxy. Such a compo- 
nent would contribute a higher fraction of the to- 
tal quasar light in the faint phases, making quasars 
appear redder than in the bright phases. Any rea- 
sonable host galaxy spectral energy distribution 
and luminosity would contribute very little light 
to the bluest bands, and would quickly be red- 
shifted out of the other bands. By a redshift of 
0.5, there should be almost no significant contam- 
ination from the host galaxies in any of the pass- 
bands. A host galaxy component cannot account 
for the wide range of redshifts over which the color 
difference is significant. 

6. Variability of Radio, X-ray, and Broad 
Absorption Line Quasars 

There is evidence from previous studies that 
the variability amplitude of quasars varies among 
different subclasses, such as those with radio 
emission (Helfand et al. 2001; Eggers, Shaffer, & 
Weistrop 2000; Garcia, Sodre, Jablonski, & Ter- 
levich 1999; Pica & Smith 1983) or broad absorp- 
tion line systems (Sirola et al. 1998). The entire 
class of highly variable blazars, for example, is 
defined in part by X-ray and radio emission (e.g. 
Ulrich, Maraschi, & Urry 1997). Here we examine 
the variability of broad subclasses of quasars in 
comparison to carefully selected control samples. 

6.1. Radio Detected Quasars 

Some SDSS quasar candidates are selected as 
optical matches to radio sources in the FIRST sur- 
vey (White, Becker, Helfand, & Gregg 1997; Ivezic 
et al. 2002). About 10% of the verified quasars in 
the sample have counterparts in the FIRST sur- 
vey. In the areas of the SDSS sample covered by 
the FIRST survey at the time the quasar candi- 
dates were selected, there are 1553 verified quasars 
with FIRST catalog matches. To test whether 
quasars with detected radio emission are more or 
less variable than those without, we have extracted 
a control sample of quasars without matches in the 
FIRST catalog. The non-radio-detected control 
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sample was constructed to have the same redshift, 
luminosity, and time lag distribution as the radio- 
detected sample, by matching each radio quasar 
with a non-radio quasar with nearly the same red- 
shift, magnitude, and time lag. The standard 
deviations of the redshift, magnitude, and time 
lag differences are er(Az) = 0.03, <r(Am) = 0.04 
mag., and ct(A(At)) = 3.6 days, and in no case 
were the differences allowed to be greater than 
0.5, 0.5 mag., and 40 days respectively. Of the 
FIRST matched quasars, 1376, 1389, and 1388 
were able to be matched with counterparts in the 
non-FIRST sample in the <?, r, and i bands respec- 
tively. Kolmogorov-Smirnov tests comparing the 
redshift, magnitude, and time lag distributions of 
the radio and control samples show that they are 
statistically indistinguishable. This also guaran- 
tees that the wavelength coverages of the samples 
are nearly identical. 

A radio-loud (not simply radio-detected) sub- 
sample and its corresponding non-radio control 
sample were also generated. Radio loudness is 
defined here as the ratio of the rest frame 5GHz 
to 4500A flux-densities (e.g. Sramek & Weedman 
1980), and a quasar is deemed radio-loud if the 
ratio is at least 100. The sample sizes of the radio- 
loud quasars and the matched radio-quiet control 
quasars are 492, 530, and 546 objects for the g, r, 
and i bands respectively. 

The matched time lag structure functions for 
the full radio and non-radio samples are shown in 
Fig. 18. There is no significant difference in the 
binned structure functions of the radio-detected 
and undetected quasars. On the other hand, the 
matched structure functions for the radio-loud 
sub-sample, shown in Fig. 19, are about 1.3 times 
higher than the non-radio sample. However, only 
the difference in the i band is significant (matched 
pair i-test probability of 1% if the samples were 
not truly different). Thus, there is marginal evi- 
dence that radio-loud quasars are more optically 
variable than radio-quiet quasars, but a larger 
sample will be needed to confirm this. 

The qualitative result that the radio-loud 
quasars arc more variable agrees with most of 
the past suggestions (Helfand et al. 2001; Eggers, 
Shaffer, & Weistrop 2000; Garcia, Sodre, Jablon- 
ski, & Terlevich 1999; Pica & Smith 1983). Since 
most blazars are radio-loud (e.g. Ulrich, Maraschi, 
& Urry 1997), the higher variability amplitude 



of radio-loud quasars may reflect a higher frac- 
tion of blazars. There is not enough information 
from this survey to reliably classify individual ob- 
jects as blazars (at the very least, more detailed 
lightcurves are needed). Further subdivision of the 
sample by finer radio loudness currently yields too 
few quasars for meaningful comparisons. In any 
case, the evidence for a correlation between radio 
loudness and UV/optical variability amplitude is 
suggestive, but not conclusive. 

6.2. X-Ray Detected Quasars 

As with radio quasars, some of the SDSS 
quasars are selected for spectroscopic follow-up 
as matches to sources in the Rosat All Sky Sur- 
vey (RASS, Voges ct al. 1999, 2000). A detailed 
analysis of RASS source matches to the SDSS 
data is given by Anderson et al. (2003). About 
5% of the verified quasars can be matched with 
RASS sources, giving about 1300 X-ray quasars 
in our sample. We constructed the X-ray and 
control samples in the same way as for the radio 
sample and its control. The numbers of matched 
objects in each of the g, r, and i bands are 1010, 
1008, and 1009 respectively. Again, Kolmogorov- 
Smirnov tests show that with respect to redshift, 
luminosity, and time lag, the two sets of samples 
are indistinguishable. 

The time lag structure functions for the X- 
ray detected and X-ray-non-detected sample are 
shown in Fig. 20. The X-ray sample is more vari- 
able than the non- X-ray control sample at time 
lags up to about 250 days, after which the differ- 
ences of variability amplitudes are much smaller. 
Overall the X-ray sample amplitudes are larger by 
a factor of 10%. The matched pair t-test prob- 
abilities of the differences occurring by chance are 
9.1%, 0.8%, and 0.4% for the g, r, and i bands re- 
spectively. The data therefore show that X-ray se- 
lected quasars are significantly more variable than 
their X-ray-non-detected counterparts. However, 
the difference may become less significant at longer 
time lags and at longer wavelengths. 

The higher X-ray variability amplitude is prob- 
ably not surprising given the high fraction of X-ray 
detected objects among blazars. As for the radio 
sample, the X-ray sample was selected purely by 
optical matches to catalog sources; no information 
about the variability of the objects in the sample 
was used beforehand. This is the first time that a 
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large X-ray selected sample of quasars has been ex- 
amined for UV/optical variability. The SDSS sam- 
ple will soon be large enough to subdivide it fur- 
ther by X-ray brightness. In the meantime, a cor- 
relation between X-ray emission and UV/optical 
variability amplitude can be claimed. We discuss 
this further in § 7. 

6.3. Broad Absorption Line Quasars 

Broad absorption line quasars (BALQSOs) are 
defined by the presence of very strong, blue-shifted 
absorption troughs in their spectra. About 5% of 
the quasars in the SDSS sample can be classified 
as BALQSOs, a fraction which is heavily redshift 
and color dependent. The largest systematically 
selected samples of BALQSOs are those of Re- 
ichard et al. (2003) and Tolca, Krolik, & Tsve- 
tanov (2002), each of which contain close to the 
same sets of objects drawn (in somewhat different 
ways) from the SDSS Early Data Release quasar 
sample. Both sets contain about 200 objects, and 
for the present purposes the samples are indistin- 
guishable; we have used the Reichard et al. (2003) 
sample since the selection process is more auto- 
mated and is likely to be used for future BALQSO 
catalogs. A control sample of non-BAL quasars 
was designed to have consistent redshift, luminos- 
ity, and time lag distributions, in the manner de- 
scribed above. The matched sample of BALQSOs 
contains 178, 189, and 190 objects in the g, r, and 
i bands respectively. 

The matched structure functions are shown in 
Fig. 21. The time lag binning is necessarily coarse 
due to the relatively small sample sizes. At the 
level of sensitivity of this sample, there is no dif- 
ference in the variability amplitudes of BAL and 
non-BAL quasars. Currently favored models of 
the BAL phenomena attribute the absorption to 
high opacity gas, either as clouds or flows, usu- 
ally viewed near the plane of an accretion disk. 
If the presence of BALs is purely a viewing an- 
gle effect, then continuum variability, if due to the 
central quasar engine is unlikely to be correlated. 
However, if variability is due to the presence of 
obscuring dust of varying attenuation crossing the 
sightline to a quasar, BAL quasars may be ex- 
pected to be more highly variable at optical and 
UV wavelengths. The issue will need to be settled 
with a larger sample, but the current results do 
not support any correlation between the presence 



of BAL features and UV/optical variability. 

7. Discussion 

To summarize our results, we have separated 
the dependence of variability on a number of pa- 
rameters, and found a power law dependence on 
time lag, anticorrelations with wavelength and lu- 
minosity, and a correlation with redshift. All of 
these relationships have been parameterized. Ra- 
dio loud and X-ray quasars also appear to be 
more variable than their quiet counterparts. There 
is currently no model of quasar continuum vari- 
ability at optical and UV wavelengths that ad- 
dresses all of the relationships described here, and 
until recently, there were virtually no quantita- 
tive predictions. Current models can be classified 
broadly into three groups: accretion disk instabil- 
ities, discrete-event or Poissonian processes, and 
gravitational microlensing. We ignore other evi- 
dence for or against the theories and describe how 
variability as an independent phenomenon may 
constrain the models. 

The Poissonian models postulate that quasar 
luminosity, or at least a significant fraction of it, 
is generated by some type of multiple discrete and 
independent energetic events, such as supernovae 
or star collisions (e.g. Terlevich, Tenorio-Tagle, 
Franco, & Melnick 1992; Torricelli-Ciamponi, 
Foellmi, Courvoisier, & Paltani 2000). The sta- 
tistical superposition of the light curves of the 
randomly occurring events determines the lumi- 
nosity at any given time. As discussed in § 4.3, the 
simplest Poissonian models predict a luminosity 
dependent power law slope of (3 = 1/2 which is 
inconsistent with our results. More detailed mod- 
els in which the event rate, energy, timescale, and 
background contribution are adjustable parame- 
ters can produce a wide range of slopes (Cid Fer- 
nandes, Sodre, & Vieira da Silva 2000), but a value 
of f3 = 1/2 is still difficult to avoid in models invok- 
ing supernovae and their remnants as the events 
(Paltani & Courvoisier 1997; Aretxaga, Cid Fer- 
nandes, & Terlevich 1997). Another apparently 
unavoidable consequence of the Poissonian models 
is that the variable luminosity component is not 
wavelength dependent (Cid Fernandes, Sodre, & 
Vieira da Silva 2000), and any color changes must 
be the result of a non- variable component (such as 
a host galaxy), which must be red to qualitatively 
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account for the wavelength correlation found here 
and in other studies. We have shown that the vari- 
able source itself is wavelength dependent, and a 
host galaxy component alone cannot account for 
the color changes. Quantitative predictions for 
the power law slope of the structure function in 
the starburst (supernova) model (Kawaguchi, Mi- 
ncshige, Umemura, & Turner 1998) range from 
7 « 0.7 — 0.9, which are also quite inconsistent 
with the value we find (7 = 0.246). Thus, based 
on predictions for the time lag, wavelength, and 
luminosity dependence of variability, current Pois- 
sonian models are inconsistent with the observa- 
tional results of this study. It remains to be seen 
if non-Poissonian processes, for example in which 
the events are not independent or random, can 
account for the observations. 

The idea that the motions of intervening matter 
along the geodesies to quasars may cause flux vari- 
ations (microlensing) was discussed as early as the 
late 1970s (Chang & Refsdal 1979), but few quan- 
titative predictions have been worked out with re- 
gard to typical quasar variability. Using the simu- 
lated microlensing light curves of Lewis, Miralda- 
Escude, Richardson, & Wambsganss (1993) and 
Schneider & Weiss (1987), Hawkins (2002) gen- 
erated structure functions which have a power 
law form with slopes ranging from 7 « 0.23 to 
7 w 0.31 (Hawkins 2002), which is consistent 
with what we find. However, the slopes depend 
upon the unknown lens mass distribution func- 
tion, velocity distribution, and source size (Wyithc 
& Turner 2001), so that a wide range of val- 
ues is possible. There is little doubt that mi- 
crolensing of quasar images does happen, and it 
has likely been detected in at least two cases, 
Q2237+0305 (Schmidt et al. 2002, and references 
therein) and Q0957+561 (Refsdal, Stabell, Pelt, & 
Schild 2000). However in each case, the quasar is 
macrolensed by a foreground galaxy, which means 
that the geodesies selectively pass through re- 
gions of relatively high density. Wyithe & Turner 
(2002) showed that the probability of microlens- 
ing by stars among single image (not macrolensed) 
sources is very small; dark matter composed of 
compact objects can improve the probability, but 
at most only about 10% of sources are expected 
to be microlensed at any given time. In addition, 
since unresolved macrolensed quasars - which ap- 
pear more luminous than they really are - are 



more likely to be microlensed, the anticorrelation 
of variability amplitude with luminosity is oppo- 
site to the trend that would be expected from mi- 
crolensing. While the deflection of light by gravity 
is achromatic, the wavelength dependence of vari- 
ability is not necessarily evidence against the mi- 
crolensing hypothesis, as long as quasar emitting 
regions are smaller and brighter with decreasing 
wavelength. It is difficult, however, to see why 
there would be any dependence on the radio or 
X-ray properties of the quasar sample. The am- 
plitude of quasar variability changes with redshift, 
but it is nearly as strong at low redshifts as it is at 
high redshifts, which also shows that microlensing 
can not be the primary cause of variability, since 
microlensing events should be extremely rare at 
low redshift. Finally, reverberation mapping stud- 
ies (Peterson 2001) show that the quasar broad 
line region varies in response to changes in the con- 
tinuum luminosity, showing that a large fraction 
of variability must be intrinsic to quasars. If some 
quasar variability is due to microlensing, it will be 
important to isolate it from the other sources since 
it has the potential to constrain the components 
of dark matter and the structure of quasars. 

It is widely accepted that quasar luminosity is 
generated through some set of processes related 
to the accretion of matter from a disk onto a 
supermassive black hole (e.g. Rees 1984). It is 
therefore natural to consider mechanisms associ- 
ated with changes in these processes as a source 
of quasar variability. Qualitatively, most schemes 
would tend to follow the trends we find here, in 
particular the anticorrelations of variability with 
wavelength and luminosity. For example, the 
disk emission spectrum of the standard optically 
thick geometrically thin accretion disk model (e.g. 
Shakura & Sunyaev 1973) is more luminous and 
bluer when the accretion rate is higher, and the 
relative luminosity change would be lower in more 
luminous objects for a given change in accretion 
rate. However, it is not known how the accre- 
tion rate would change nor how the resultant lumi- 
nosity changes would propagate through the disk. 
While there has been much theoretical discussion 
of the various possibilities for the emission mecha- 
nisms and their instabilities (e.g. Wallinder, Kato, 
& Abramowicz 1992; Schramkowski & Torkels- 
son 1996), there are currently few quantitative 
predictions which can be compared to observa- 
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tional results. Kawaguchi, Mineshige, Umemura, 
& Turner (1998) generated structure functions for 
the cellular-automaton model for disk instability, 
and found power law forms with slopes ranging 
from 7 = 0.41 to 0.49. While this range is incon- 
sistent with our results, the model is necessarily 
simplified and a number of assumptions need to 
be made. The complexity of possible accretion 
disk (or jet) instability models is likely what has 
prevented more quantitative predictions. Disk in- 
stability models are clearly promising, but as yet 
it is difficult to compare them to the observations. 

In summary, the weight of the observational ev- 
idence seems to disfavor gravitational microlensing 
and generic Poissonian processes as the primary 
source of quasar variability. Accretion disk insta- 
bility models have yet to be adequately developed 
quantitatively for direct comparison with our re- 
sults. It is also plausible that a combination of 
sources produce variations in quasar lightcurves, 
and no single model can be completely eliminated 
at this time. 

8. Conclusions 

We have examined the ensemble broadband 
photometric variability of a very large and homo- 
geneous sample of quasars from the SDSS - the 
largest sample ever used to study variability. The 
three-band spectrophotometry of each object was 
compared directly to the imaging photometry ob- 
tained at an earlier epoch. Because of the large 
number of objects and wide coverage of param- 
eter space, the dependences of variability ampli- 
tude on time lag, luminosity, wavelength, and red- 
shift were able to be disentangled for the first time. 
The variability amplitude increases with time lag 
(up to about two years) as a power-law with a 
slope of 7 = 0.25. In terms of the variability 
amplitude, more luminous quasars are less vari- 
able, shorter wavelengths arc more variable, and 
more distant quasars are somewhat more variable; 
all of these relationships are parameterized. Ra- 
dio loud quasars appear to be more variable than 
their radio quiet counterparts, and quasars with 
detectable X-ray emission (in the ROSAT survey) 
are more variable than those without. It is diffi- 
cult to explain the results in the context of models 
involving discrete events (Poissonian models) and 
gravitational microlensing. Accretion disk insta- 



bility models are promising, but more quantita- 
tive predictions are needed to test them against 
the observational results. 
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Fig. 1. — The uncorrected spectroscopic minus photometric magnitudes vs. imaging PSF magnitudes for 
stars observed on the same spectral plates as the quasars (left). Results for all three bandpasses are shown. 
The curves show the binned median trends and the upper and lower 68.3% confidence envelopes. The right 
side panels show the uncorrected magnitude difference histograms. The 68.3% confidence half-widths, a, are 
given for a spectroscopic S/N of 10. 
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Fig. 2. — The corrected spectroscopic minus photometric magnitudes vs. spectral signal-to-noise ratio 
for stars observed on the same spectral plates as the quasars (left). Results for all three bandpasses are 
shown. The curves are the fits to the 68.3% confidence half-width envelopes. The right side panels show 
the magnitude difference histograms for each band. The 68.3% confidence half-widths, a are given for a 
spectroscopic S/N of 10. 
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Fig. 3. — The corrected spectroscopic minus photometric magnitudes vs. spectral signal-to-noise ratio for 
the quasars (left). Results for all three bandpasses are shown. The curves are the binned 68.3% confidence 
half-width envelopes. The right side panels show the magnitude difference histograms for each band. The 
68.3% confidence half-widths, a, are given for a spectroscopic S/N of 10. 
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Fig. 4. — Magnitude difference (uncorrected for measurement uncertainties) vs. rest frame time delay in each 
of the three photometric pass bands. The binned points show the mean values while the error bars show the 
root-mean-square deviations divided by the square root of the number of objects in a bin. 
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Fig. 5. — Magnitude difference (uncorrected for measurement uncertainties) vs. absolute i band magnitude 
in each of the three photometric pass bands. The binned points show the mean values while the error bars 
show the root-mean-square deviations divided by the square root of the number of objects in a bin. 
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Fig. 6. — Magnitude difference (uncorrected for measurement uncertainties) vs. rest wavelength in each of 
the three photometric pass bands. The binned points show the mean values while the error bars show the 
root-mean-square deviations divided by the square root of the number of objects in a bin. 
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Fig. 7. — Magnitude difference (uncorrected for measurement uncertainties) vs. rcdshift in each of the three 
photometric pass bands. The binned points show the mean values while the error bars show the root-mean- 
square deviations divided by the square root of the number of objects in a bin. 
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Fig. 8. — Quasar structure functions for each of the three pass bands, color coded by band. No accounting 
has been made for any other variability dependencies, such as luminosity, wavelength, or redshift. Single 
power law fits to the data are also shown. 
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Fig. 9. — Projected parameter values for all of the quasars. This effectively shows the selection function in 
the parameter space given by rest frame time lag, At, rcdshift, z, absolute i-band magnitude, Mi, and rest 
wavelength, A reS i . The passbands are indicated when necessary. For clarity, only the r band data are shown 
for the middle two plots. 
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Fig. 10. — Variability as a function of rest frame time lag for independent slices of data (upper left, upper 
right, and lower left panels). Colors indicate redshift, with redder colors showing data from higher redshift 
ranges, according to the color key. Solid and open points show results from slices with higher and lower 
luminosities respectively, within the same redshift range. Results from the three photometric bands are given 
separately, which effectively restricts the rest wavelength ranges to small values in each slice. The lower right 
panel shows all of the scaled data points, along with the best fit power law for variability vs. At. 
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Fig. 11. — Variability as a function of absolute magnitude for independent slices of data (upper left, upper 
right, and lower left panels). Colors indicate redshift, with redder colors showing data from higher redshift 
ranges, according to the color key. Solid and open points show results from slices with longer and shorter 
time lags respectively, within the same redshift range. Results from the three photometric bands are given 
separately, which effectively restricts the rest wavelength ranges to small values in each slice. The lower right 
panel shows all of the scaled data points, along with a best fit generalized Poissonian function (solid) and 
Poissonian function with the power law index fixed at 1/2 (dashed). 
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Fig. 12. — Variability as a function of rest frame wavelength for independent slices of data. Colors indicate 
average rest frame time lag, with redder colors showing longer time lags, according to the color key. Solid 
and open points show results from slices with more luminous and less luminous quasars respectively, within 
the same redshift and time lag ranges. Results from the six redshift slices are shown in separate panels. 
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Fig. 13. — Scaled variability amplitude points from Fig. 12 as a function of rest frame wavelength. The line 
shows the best fit exponential function as described in the text. 
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Fig. 14. — The triangle shows the region ol the wavelength-absolute magnitude plane that is simultaneously 
well covered by quasar observations in all three photometric bands. The vertical lines - separated by 100, 
200, 400, or 800A - show the bin sides used to separate the quasar sample by wavelength. 
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Fig. 15. — Variability as a function of redshift for independent sets of data. Sets of three connected points 
show the variability amplitude for data sets with similar distributions of rest wavelength, absolute magnitude, 
and time lag, but different average redshifts. Colors correspond to the average absolute magnitude for each 
set, with bluer colors indicating brighter absolute magnitudes, according to the color key. The point sizes 
correspond to average time lags, according to the point size scale. 
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Fig. 16. — Scaled variability amplitude as a lunction of redshift. The best fit line (to linear redshift and 
variability amplitude) is also shown. 
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Fig. 17. — Average faint phase vs. bright phase colors of quasars as a function of redshift (left). The 
structure in the plots is mainly due to spectral features redshifting into and out of the various passbands. 
The plots on the right side show the color difference between the bright and faint phases as a function of 
redshift, along with the uncertainties in each bin. A more negative value implies a bluer spectrum in the 




Fig. 18. — Radio-detected (solid lines and points) vs. undetected (dashed lines and open points) structure 
functions. 



37 




60 80 100 200 

Quosor Rest Frome Time Log, At (doys) 



Fig. 19. — Radio- loud (solid lines and points) vs. radio-quiet (dashed lines and open points) structure 
functions. 
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Fig. 20. — X-ray-detected (solid lines and points) vs. undetected (dashed lines and open points) structure 
functions. 
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Fig. 21. — BALQSO (solid lines and points) vs. non-BALQSO (dashed lines and open points) structure 
functions. 
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Table 1 

Parameter values for fits to the binned structure functions. 



Power Law 


Band 


At (days) 7 


x 2 


9 


9.90 ± 6.49 x 10 4 0.293 ±0.030 


1.51 


r 


7.05 ± 4.23 x 10 4 0.336 ±0.033 


1.20 


i 


1.66 ± 1.32 x 10 s 0.303 ±0.035 


4.20 


Exponential 


Band 


V (mag) At (days) 


x 2 


9 


0.168 ± 0.005 51.9 ±6.0 


20.5 


r 


0.155 ±0.006 74.7 ±8.9 


24.8 


i 


0.139 ± 0.005 62.6 ±8.3 


39.3 



Table 2 

Parameter values for fits to the scaled variability time lag dependence. 



Power Law 


Atq (days) 7 


x 2 


5.36 ± 1.46 x 10 5 0.246 ±0.008 


299.0 


Exponential 


Vo At (days) 


x 2 


0.144 ±0.001 40.4 ±1.4 


443.4 
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